The star-formation history of the universe - an infrared 

perspective 
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ABSTRACT 

A simple and versatile parameterized approach to the star formation history allows a quanti- 
tative investigation of the constraints from far infrared and submillimetre counts and background 
intensity measurements. 

The models include four spectral components: infrared cirrus (emission from interstellar dust), 
an M82-like starburst, an Arp220-like starburst and an AGN dust torus. The 60 luminosity 
function is determined for each chosen rate of evolution using the PSCz redshift data for 15000 
galaxies. The proportions of each spectral type as a function of 60 fim luminosity are chosen for 
consistency with IRAS and SCUBA colour-luminosity relations, and with the fraction of AGN as 
a function of luminosity found in 12 fim samples. The luminosity function for each component 
at any wavelength can then be calculated from the assumed spectral energy distributions. With 
assumptions about the optical seds corresponding to each component and, for the AGN com- 
ponent, an assumed dependence of the dust covering factor on luminosity, the optical and near 
infrared counts can be accurately modelled. High and low mass stars are treated separately, since 
the former will trace the rate of star formation, while the latter trace the cumulative integral of 
the star formation rate. 

A good fit to the observed counts at 0.44, 2.2, 15, 60, 90, 175 and 850 /zm can be found with 
pure luminosity evolution in all 3 cosmological models investigated: f2 D = 1, f2 D = 0.3 (A = 0), 
and fl = 0.3, A = 0.7. 

All 3 models also give an acceptable fit to the integrated background spectrum. Selected 
predictions of the models, for example redshift distributions for each component at selected 
wavelengths and fluxes, are shown. The efect of including an element of density evolution is also 
investigated. 

The total mass-density of stars generated is consistent with that observed, in all 3 cosmological 
models. 

Subject headings: infrared: galaxies - galaxies: evolution - star:formation - galaxies: starburst - cosmol- 
ogy: observations 



1. Introduction 

In this paper I investigate what source counts 
and the background radiation at infrared and 
submm wavelengths can tell us about the star for- 
mation history of the universe. Madau et al (1996) 
showed how the ultraviolet surveys of Lilly et al 
(1996) could be combined with information on uv 
dropout galaxies in the Hubble Deep Field to give 
an estimate of the star formation history from z 



= 0-4 (see also Madau et al 1998). However these 
studies ignored what is already well-known from 
far infared wavelengths, that dust plays a major 
role in redistributing the energy radiated at visible 
and uv wavelengths to the far infrared. 

Subsequent to the Madau et al analysis several 
groups of authors have argued that the role of dust 
is crucial in estimates of the star formation rates 
at high redshifts. Rowan-Robinson et al (1997) 
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derived a surprisingly high rate of star formation 
at z — 0.5-1 from an ISO survey of the Hubble 
Deep Field (HDF). Subsequent ISO estimates by 
Flores et al (1999) appear to confirm the need to 
correct for the effects of dust in estimates of star- 
formation rates. 

Large extinction correction factors (5-10) were 
derived at z = 2-5 by Meurer et al (1997, 1999). 
Pettini et al (1998), Dickinson (1998) and Steidel 
et al (1999) found that estimates of star formation 
rates derived from HDF and other Lyman drop- 
out studies do indeed need to be corrected upwards 
by a substantial factor (2-7). Thus the high star 
formation rates deduced from the ISO-HDF data 
(Rowan-Robinson et al 1997) are not as out of line 
as at first appeared. 

Many authors have attempted to model the 
star formation history of the universe ab initio (eg 
Kauffmann et al 1993, Cole et al 1994, Pei and 
Fall 1995, Baugh et al 1998, Blain et al 1999c, 
Somerville et al 2000, Kauffmann and Haehnclt 
2000). There is an even longer history of attempts 
to find empirical models for the evolution of the 
starburst galaxy population (eg Franceschini et al 
1991, 1994, 1997, Blain and Longair 1993, Pear- 
son and Rowan-Robinson 1996, Guiderdoni et al 
1997, 1998, Dwek et al 1998, Blain et al 1999a, 
Dole et al 2000). The study of Blain et al (1999a) 
attempts a more systematic approach with a pa- 
rameterized set of models but still treat the spec- 
tral energy distributions (seds) of the galaxies in 
a simplistic manner. In the present paper I report 
the results of a parameterized approach to the star 
formation history of the universe, which allows a 
large category of possible histories to be explored 
and quantified. The parametrized models can be 
compared with a wide range of source-count and 
background data at far infrared and submillimtre 
wavelengths to narrow down the parameter space 
that the star formation history can occupy. The 
approach is similar to that of Guiderdoni et al 
(1998) and Blain et al (1999a), but differs in key 
respects outlined below. 

An important development for understanding 
the role of dust in the star formation history of 
the universe has been the detection of the submm 
background at wavelengths from 140-850 /im using 
COBE data (Puget et al 1996, Schlegel et al 1998, 
Hauser et al 1998, Fixsen et al 1998, Lagache et al 
1999). There is still disagreement about the actual 



value of the background intensity by a factor of 2 
at 850 /im and 2.5 at 140 /um. 

Franceschini et al (1997), Guiderdoni et al 
(1998), Blain et al (1999a), Dwek et al (1998) 
and Gispert et al (2000) have presented models to 
account for this background radiation. The mod- 
els of Guiderdoni et al are able to account for the 
350 and 850 (im background, but fail to fit the 
observed background at 140 and 240 /zm, even if 
the lower values of Fixsen et al are adopted. The 
models of Franceschini et al (1997) invoke a new 
population of heavily obscured high redshift star- 
bursts, designed to account for the formation of 
spiral bulges and ellipticals. In this paper I test 
whether the submillimetre background and counts 
can be understood in terms of a single popula- 
tion of evolving star-forming galaxies. Blain et 
al (1999a) put forward a range of models which 
they argue are consistent with all available counts 
and background radiation. I shall argue below 
that, amongst other problems, their most success- 
ful models involve an unphysical discontinuity in 
the star formation rate. Of the two most suc- 
cessful models of Dwek et al (1998), one ('RR') is 
taken from an earlier version of the present work, 
and the other ('ED') involves a spike in the star 
formation history. 

The layout of the paper is as follows. In section 
2, I introduce a parameterized characterization of 
the star formation history which is both simple 
and versatile. Section 3 discusses the 60 lumi- 
nosity function and evolution rate as derived from 
the IRAS PSCz galaxy redshift survey. Section 
4 discusses the assumed spectral energy distribu- 
tions of galaxies as a function of luminosity, used 
to predict the luminosity functions at other wave- 
lengths. Section 5 gives combined least squares 
fits to the counts at 60 and 850 /an, and the back- 
ground intensity at 140-750 /zm in 3 selected cos- 
mological models. Section 6 discusses the prop- 
erties of the best-fitting models at a wide range 
of wavelengths. The effect of including some den- 
sity evolution in an tt D = 1 model is discussed in 
section 7. Finally section 8 gives discussion and 
conclusions. 

A Hubble constant of 100 km/s/Mpc is used 
throughout. 
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2. Parametrized approach to star forma- 
tion history 

In this paper I present a parameterized ap- 
proach to the problem, investigating a wide range 
of possible star formation histories for consistency 
with counts and background data from the ultra- 
violet to the submm. 

The constraints we have on the star formation 
rate, 0* (t) are that 

(i) it is zero for t = 

(ii) it is finite at t = t 

(iii) it increases with z out to at least z = 1 ( 
and from (i) must eventually decrease at high z) . 

A simple mathematical form consistent with 
these constraints is 

MtVMto) = cxp Q(i-(W) {t/t ) p (i) 

where P and Q are parameters 

( P > to satisfy (i), Q > to satisfy (iii)). I 
assume that 4>*(t) = for z > 10. 

Equation (1) provides a simple but versatile pa- 
rameterization of the star formation history, capa- 
ble of reproducing most physically realistic, single- 
population scenarios. Figure 1 shows models of 
type eqn (1) (a) chosen to mimic the Pei and Fall 
(1995) models based on analysis of the history of 
heavy element formation, and also the Cole et al 
(1994) semianalytic simulations of the star forma- 
tion history in a cold dark matter (CDM) scenario, 
(b) chosen to be consistent with the strong evolu- 
tion rate seen by Lilly et al (1996) in the CFRS 
survey, (c) fitted to far infrared and submillimetre 
source-counts and background intensity in differ- 
ent cosmological models (see below). Models of 
the form (1) will not, however, be able to repro- 
duce the very sharply peaked scenario ('ED') of 
Dwek et al (1998), or the two-population model of 
Franceschini et al (1997). In the scenarios mod- 
elled by eqn (1) the formation of bulges and el- 
lipticals are part of an evolving spectrum of star- 
formation rates. 

The physical meaning of the parameters is as 
follows. The exponential term describes an expo- 
nential decay in the star formation rate on a time- 
scale t a /Q, which can be interpreted as a prod- 
uct of the process of exhaustion of available gas 
for star formation (a competition between forma- 



tion into stars and return of gas to the interstellar 
medium from stars) and of the declining rate of 
galaxy interactions and mergers at later epochs. 
This parameter is the same as that used in the 
galaxy sed models of Bruzual and Chariot (1993). 
The power-law term represents the build-up of the 
rate of star formation due to mergers of fragments 
which are incorporated into galaxies. P measures 
how steeply this process occurs with time. The ra- 
tio P/Q determines the location of the peak in the 
star formation rate, t pea k, since t pea k/t = P/Q. 

A very important assumption in the present 
work is that the star formation rate should vary 
smoothly with epoch. Several earlier assumptions 
have assumed, for mathematical convenience, dis- 
continuous changes of slope in the star formation 
rate (eg Pearson and Rowan-Robinson 1996, the 
Blain et al 1999a 'anvil' models). Such discontinu- 
ities are highly unphysical and I have eliminated 
them in this work. The assumption above that 
(j)*(t) = for z > 10 represents just such an un- 
physical discontinuity in eqn (1) and I have inves- 
tigated the effect of replacing the power-law term 
in eqn (1) by 

(t/*o) P - (ti/to) P 

where tj is the epoch at which star forma- 
tion commences. This would then constitute a 
third parameter in the models, but one which has 
very little effect on the results, provided P > 
l,z(ti) » 1. The Blain et al (1999a) 'peak-4' 
and 'peak-5' models assume that there is no star- 
formation at z < 4, 5, respectively, and this consti- 
tutes a strong discontinuity within the observable 
range of redshifts, so would again not be consis- 
tent with the philosophy adopted here. The Blain 
et al (1999a) 'peak-10' model is reasonable well 
approximated by (P,Q) = (1.9, 10.15). 

We might expect that the cosmological model 
could have a significant effect on the relationship 
between predicted counts and predicted back- 
ground intensity, since the latter is sensitive 
to how the volume element and look-back time 
change with redshift. 

To test this I have explored models with (a) A = 
0, for which all the required formulae are analytic 
(specifically, here, models with Q = 1, 0.3), and 
(b) k = 0, for which some of them are (specifically, 
A = 0.7, fl a = 0.3). Relevant formula are given in 
Appendix A. 



3 



3. 60 /im luminosity function and evolu- 
tion derived from IRAS PSCz survey 

Given an assumed (P,Q) I then determine the 
60 fim luminosity function, using the IRAS PSCz 
sample (Saunders et al 1999). I fit this with the 
form assumed by Saunders et al (1990) 

7](L) = C4L/U.) 1 - a e- ^ t °9^ 1 + L / L '1/^ 2 (2) 

and find that the luminosity function parame- 
ters can be well approximated as follows: 

for fl = 1: 

a = 0.744-0.033 W 

log w L* = 8.50 + 0.05 W - 2 log w (H o /100) 
where W = Q - 1.125 P 
for Q = 0.3(A = 0): 
a = 0.748-0.028 W 

log w L* = 8.46 + 0.05 W - 2 log w (H o /100) 
where W = Q - 1.10 P 
for A = 0.7(fc = 0): 
a = 0.783-0.030 W 

log w L* = 8.39 + 0.05 W - 2 log w (H o /100) 
where W = Q - 1.067 P. 

I have assumed fixed values for a = 1.09, C* = 
0.027 (if o /100) 3 . It is not clear that any previous 
studies have correctly taken account of the need to 
change the 60 /im luminosity function as the rate 
of evolution is varied. The study of Guiderdoni et 
al (1998) explicitly violates the known constraints 
on the 60 fim luminosity function at the high lu- 
minosity end and as a result the models predict far 
too many high rcdshift galaxies at a flux-limit of 
0.2 Jy at 60 /jm, where substantial redshift surveys 
have already taken place. Blain et al (1999a) state 
that they have determined the luminosity function 
from consistency with the 60 /jm counts, but this 
process does not automatically give detailed con- 
sistency with existing redshift surveys. 

We can also use the PSCz data to determine 
a range of consistency for (P,Q), using the V/V m 
test. The predicted uncertainty in < V/V m > for 
a population of n galaxies, (12n) -1 / 2 can be used 
to assign a goodness of fit for each set of (P,Q) 



values. Fig 7 shows the locus of the best-fitting 
models in the P-Q plane (for which < V/V m > = 
0.5), together with the ±1 — a region, for Q D = 1. 

4. Assumed infrared and submm spectral 
energy distributions 

To transform this 60 /im luminosity function to 
other wavelengths, we have to make an assumption 
about the spectral energy distributions. I have ex- 
plored a variety of assumptions about these seds: 
(a) a single sed at all luminosities representative of 
starbursts; (b) composite seds which are a mixture 
of four components, a starburst component, a 'cir- 
rus' component, an Arp 220-like starburst, and an 
AGN dust torus (cf Rowan-Robinson and Craw- 
ford 1989, Rowan-Robinson and Eftstathiou 1993, 
Rowan- Robinson 1995), with the proportions of 
each depending on 60 /jm luminosity (Rowan- 
Robinson 1999a, b). Neither of these approaches 
gave satisfactory results and it was not possible to 
find a simultaneous fit to all the far infrared and 
submillimtre counts and background spectrum in 
any cosmological model. Models with luminosity- 
dependent composite seds could be found which 
fit the far infrared and submillimetre counts and 
850 (im background, but they failed to account for 
the 140-350 /jm background. Finally I have de- 
rived counts and background spectrum separately 
for each of the four components and then summed. 

Thus 



4 r rz(L 60 ,S v M) 

(L 60 )ti(L 60 )dL 60 / (dV/dz)dz 
i=i J J ° 

(3) 

where lgL v = lgL 6Q + lg[L„/L 60 } 

= 6.18+lgS u +lg(u/u eQ )+2lgDi um (z)+k(z, v) + 
lg[Mz)/M0)l 

L Vl Lqq are in solar units, S v is in Jy, Di um is 
in Jy, and the k-correction 

k(z,v)=lg[L v (v(l + z))/L v (u)]. 
For small z, 

dN v (S„) 
7760 (Leo) (WS,) 3 / 2 

^ U{Lqo)(L v / L m ) z J 2 
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(4). 

This approach finally gave satisfactory fits to 
all available data. It also allows a correct deter- 
mination of the redshift distribution of each type 
separately as a function of wavelength and flux- 
density, and the proportion of each type contri- 
bution to the counts at any flux-density or to the 
background. 

I have used the latest predictions for infrared 
seds of these components by Efstathiou et al 
(2000), and added near ir-optical-uv seds corre- 
sponding to an Sab galaxy (Yoshii and Takahara 
1988) for the cirrus and an HII galaxy (Tol 1924- 
416, Calzetti and Kinney 1992) for the starburst 
component, respectively. The starburst model is 
a good fit to multiwavelength data for M82 and 
NGC6092 (Efstathiou et al 2000), and also to far 
infrared and submillimctre data for luminous star- 
bursts (Rigopoulou et al 1996). 

The normalization between far infrared and 
optical-uv components is determined by L(60/zm)/ 
L(0.8/zm) = 0.15 for the cirrus component, 5.3 
for the starburst component, and 49 for the Arp 
200 component. For the AGN component I as- 
sume that L(10/im)/L(0.1(/im) = 0.3 (cf Rowan- 
Robinson 1995) for the most luminous AGN, and 
that this ratio increases with decreasing luminos- 
ity to account for the fact that the mean cover- 
ing factor is higher at lower luminosities. So lg 
L(10/mi)/L(0.1(/mi) = -0.52 + 0.1*(14.0-log L60). 

For the cirrus component I have, somewhat 
arbitrarily, divided the optical sed into a contri- 
bution of young, high-mass stars (dominating at 
A < 0.3/tm) and a contribution of old, low-mass 
stars (dominating at A > 0.3/zm) (see Fig 5). The 
former are assumed to trace the star formation 
rate, and so participate in the strong evolution of 
the form cqn (1), but the latter trace the cumula- 
tive star formation up to the epoch observed, ef- 
fectively a negative evolution with increasing red- 
shift. This treatment, though approximate, al- 
lows a reasonable prediction of the K- and B- 
band counts. The two components in Fig 5 can 
be modelled, assuming L cx M 3 , blackbody seds 
with T cx L 1 / 2 , and with a Salpeter mass func- 
tion, with mass-range 0.1-1 M© for the low-mass 
star component, and 8-40 M for the high-mass 
star component. 

The proportions of the four components (at 60 



/im) clS Si function of luminosity (table 1) have 
been chosen to give the correct mean relations 
in the S(25)/S12), S(60)/S(25), S(100)/S(60) and 
S(60)/S(850) versus L(60) diagrams (Figs 2-4). 
Where predictions are being compared with IRAS 
12 /zni data, or (later) with ISO 15 /im or 6.7 tt 
counts, account is taken of the width of the rele- 
vant observation bands by filtering with a top-hat 
filter of appropriate half-width. (0.23, 0.26 and 
0.16 at 6.7, 12, and 15 /im respectively). Other- 
wise observations were assumed to be monochro- 
matic. The relative proportion of the 60 /im emis- 
sion due to AGN dust tori as a function of L(60) 
is derived from the luminosity functions given by 
Rush et al (1993) (see Fig 10). The resulting 
mean seds as a function of L(60) are shown in Fig 
6. Luminosity functions at different wavelengths 
and the infrared bolometric luminosity functions 
are shown in Figs 8-12. There is good agree- 
ment with measured luminosity functions at wave- 
lengths from 0.44-850 /im. The agreement with 
the luminosity function derived at 60 /im from 
PSCz data is unsurprising, since the luminosity 
function parameters were determined from these 
data. The fit to the 850 /im luminsity function of 
Dunne et al (2000) is more impressive, since the 
transformation from 60 to 850 /im is based only on 
choosing the U(L§q) to give the correct ridge-line 
in Fig 4. The fit to the non-Seyfert data of Rush 
et al (1993) at 12 /im is also impressive, again aris- 
ing only from the choice of ti(L§o) to be consistent 
with the colour-luminosity diagrams in Fig 3. The 
fit to the Seyfert luminosity function at 12 /im is 
not fortuitous, since the ij(L 60 ) for the AGN com- 
ponent were chosen to give this agreement. It is 
impressive that luminosity functions derived in the 
far infrared can fit the data at 0.44 /im (B band): 
the only freedom in the models to fit the B-band 
luminosity functions and counts is the amplitude 
of the optical sed relative to the far infrared. 

Clearly it will be important to have submil- 
limctre data for a wide range of normal and ac- 
tive galaxies to test and improve these seds. But 
the approach of using accurate radiative trans- 
fer models, with realistic assumptions about dust 
grains, which have been verified with observations 
of known galaxies, seems superior to modelling 
the sed as a blackbody with power-law dust grain 
opacity in which the dust temperature is treated 
as a free parameter (as in Blain et al 1999a). The 
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latter approximation can only be valid for rest- 
frame wavelengths greater than 60 /im, ie for ac- 
curate prediction of counts and background inten- 
sities at wavelengths > 200 /im. Useful predictions 
can certainly not be made at 15 /im without ex- 
plicit treatment of PAHs. These criticisms do not 
apply to the studies of Guidcrdoni et al (1998), 
Dwek et al (1998), Xu ct al (1998), Dole et al 
(2000), whose assumed seds are similar to those 
used here (but do not include the detailed depen- 
dence on luminosity required for consistency with 
Figs 2-4). 

I have also assumed that the same luminosity 
evolution function should be applied to the whole 
60 /im function, ie to AGN, 'cirrus' and 'starburst' 
components. I have investigated the effect of mak- 
ing the switchover in proportions of different types 
of component at a fixed luminosity, so that there 
is in effect a strong increase in the proportion of 
galaxies that are starbursts (or contain AGN) with 
redshift. However this did not permit a fit to all 
the available data, so I have assumed that the Ta- 
ble 1 proportions relate to the luminosity at zero 
redshift and that all these luminosities are sub- 
ject to the evolution. This is an approach which 
works, but models with more complex evolution 
than used here are clearly possible. 

A substantial part of the illumination of the 
cirrus component in spiral galaxies is due to re- 
cently formed massive stars, part of whose light 
escapes directly from star-forming regions despite 
the high average optical depth in these regions. In 
the starburst models of Efstathiou et al (2000), 
this corresponds to the late stages of their star- 
burst models. If the typical starburst luminosity 
was greater in the past then the emission from 
interstellar dust in the galaxy would also be corre- 
spondingly greater. However the application of a 
uniform evolution rate to both the starburst and 
cirrus components of the far infrared emission is 
not strictly self-consistent. At low redshifts the 
cirrus is illuminated by both low-mass and high- 
mass stars, with both contributing comparable en- 
ergy inputs. At early epochs the illumination will 
be predominantly by high-mass stars, since there 
would have been far fewer low-mass stars in the 
galaxy. So the evolution of the cirrus component 
will be complex, depending on the star-formation 
history, the evolution of metallicity and the pos- 
sibly varying dust geometry in the galaxy. The 



evolution of the seds of each component with red- 
shift, particularly of the cirrus component, may be 
an important factor in understanding the counts 
and background, but is not treated here (see dis- 
cussion in section 8). 

It is possible that the evolution of AGN differs 
from that of starburst at z > 3, but this will have 
little effect on the far infrared counts and back- 
ground (there could be a significant effect at 15 
^m, which will be worth further study). 

Elliptical galaxies are not treated explicitly, 
though their star formation rate must have been 
much greater in the past than at present. I am 
assuming that ellipticals are quiescent starburst 
galaxies, that their star formation proceeded in 
much the same way as we see in current live star- 
bursts, and that their star formation episodes are 
part of the evolution history quantified here. We 
have to think of this history as a series of short- 
lived fireworks taking place in different galaxies at 
different times. Similarly this approach does not 
track the different spiral types separately, but only 
as a global average at each epoch. 

5. Combined fits to 60, 175 and 850 /im 
counts, and 140-750 /im background, 
and determination of P,Q 

I can now predict the counts and background 
intensity at any wavelength and by comparing 
with observed values, constrain loci in the P-Q 
plane. To determine (P,Q) for any given cosmo- 
logical model, I combine the constraints found at 
60 /im from the PSCz (section 4 above) with con- 
straints from (1) deep counts at 60 /im (50 mJy), 
(2) the observed source-counts at 850 /im at 1 and 
4 mJy, (3) the background intensity at 140, 350 
and 750 /im. Figure 7 shows the constraints from 
the V/V m test and the 850 /im counts for the o = 
1 case. The other constraints are less important in 
specifying P and Q. For all 3 cosmological models, 
values of (P,Q) can be found which give a satisfac- 
tory fit to all the available data. As emphasized 
above, this outcome does depend strongly on the 
assumptions made about the seds. 

An important constraint on the models is that 
the total mass of stars produced in galaxies should 
be consistent with the mass of stars observed, 
ft* = 0.003 ± 0.0009/1- 1 (Lanzetta et al 1996), 
and that it should be less than the total density 
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of baryons in the universe, 0* < Q& = 0.0125 ± 
0.0025^ 2 (Walker et al 1991). 

We can calculate the total mass-density of stars 
from the 60 /im luminsity density using eqn (7) 
of Rowan- Robinson ct al (1997), modified to take 
account of the latest Bruzual and Chariot galaxy 
evolution models (Madau et al 1998, Rowan- 
Robinson 2000), from 

= lQ- 1113 Zh- 2 l 6Q (t o /10Gyr) (5) 

where is the luminosity density in solar lu- 
minosities per Mpc 3 , 

S = S£[Mt)/Mto)]d(t/t ), 

and the assumed fraction of opt-uv light being 
radiated in the far infrared has been assumed to 
be e = 2/3. 

Table 4 gives the values of (P,Q) which provide 
the best fit to the far infrared and submillimctrc 
counts and background, for each of the 3 cosmo- 
logical models considered, and the corresponding 
values of Iqq, £ and £1*, for an assumed age of the 
universe t a = 13 Gyr. The predicted values of 
are consistent with the oberved value in each 
case. Estimating 51* from the young stellar com- 
ponent at 2800 A or from the K-band luminos- 
ity density (with an assumed mass-to-light ratio) 
also give consistent results for an assumed Salpeter 
mass-function. 

6. Predicted counts and integrated back- 
ground spectrum from uv to submm 

Figure 13-19 shows the predicted source-counts 
in the 3 selected cosmological models, at 850, 175, 
90, 60, 15, 2.2 and 0.44 /mi. The agreement with 
observations at infrared and submillimetre wave- 
lengths is extremely impressive. Although the fits 
at 850 and 60 /mi have been ensured by the least- 
squares procedure for determining P and Q, the 
fits at 175, 90 and 15 /i are simply a consequence 
of the assumed seds and the choice of the U^Lqq). 
There is not much difference between the predic- 
tions of the 3 cosmological models at 60-850 /zm. 
At 15 /im there is a difference between the models 
in the predicted numbers of sources at fluxes be- 
low 100 /iJy. The proportion of AGN dust tori at 
12 /i agrees well that the data of Rush et al (1993) 
(15% brighter than 0.4 Jy). Fig 17 shows that the 
proportion of AGN at 15 /im is reasonably con- 
stant (15-20%) for fluxes brighter than 3 mJy, but 



is predicted to fall rapidly towards fainter fluxes. 

The fits to the counts at B and K are surpris- 
ingly good, though the model does not at present 
have ingredients capable of accounting for the very 
faint K- and B-band galaxy counts. This could 
be provided either by a measure of density evo- 
lution (see below) or by steepening the faint-end 
luminosity function at z > 1, which could in ci- 
ther case be attributed to a population of galax- 
ies that had merged into present-day galaxies. 
The B-band galaxy counts are dominated by the 
(negatively evolving) low-mass star component for 
/ 3 i S'(0.44^m) > -4.7 (Jy) (B AB < 20.7) and by 
the (positively evolving) high-mass star compo- 
nent at fainter fluxes. In the I-band (0.8 /zm) the 
crossover is at 23.7 m. The predicted median red- 
shift for Iab < 22.5 is 0.49, in good agreement 
with the results of the CFRS survey (Lilley et al 
1995), but for B < 22, the predicted median z is 
0.54, rather higher than the ^0.3 observed by Col- 
less et al (1993), suggesting that we do not have 
quite the correct luminosity function or evolution 
for the blue galaxy population. 

Fig 20-23 show redshift distributions at selected 
wavelengths and fluxes (complete output of the 
models is available at http://icstar5.ph.ac.uk/ 
~mrr/countmodels). The median redshift at 
S(850 /zm) > 2 mJy (Fig 20) is predicted to be 
2.25, significantly deeper than the prediction for 
S(0.44 fim) > 0.1 /zJy {B AB < 26.4 m.), 1.63 (Fig 
23). Only for Bab < 29 do we reach the same 
predicted median z as for 850 /im. This shows the 
power of the 850 /im surveys and also the difficulty 
there will be in identifying the sources detected. 

Figure 24 shows the predicted integrated back- 
ground spectrum for the 3 cosmological mod- 
els compared with the submillimetre observations 
summarized in Table 3, the optical and near in- 
frared estimates of Pozzetti et al (1998) and the 
15 ^m estimate of Sergeant et al (2000). Dwek 
ct al (1998a) have argued by extrapolation from 
longer wavelengths that the 100 /xm background 
should lie in the range 5-34 nWmT^sr^ 1 , but as 
this is not a true measurement of the background 
I have not used it here (my models lie within the 
estimated range). Dwek et al (1998b), Gorjian 
et al (2000) and Wright and Reese (2000) have 
given estimates of the background at 2.2 and 3.5 
/im, which appear to be a factor 2 or more higher 
than estimates derived from summing deep galaxy 
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counts (Pozzetti et al 1998). This would seem to 
imply either incomplete subtraction of foreground 
emissions or some additional background uncon- 
nected with galaxies and I have not used these 
estimates. 

The models shown in Fig 24 are consistent with 
the observations, though the predictions of the o 
= 1 model arc on the low side, while those of the 
A = 0.7 model are on the high side. Figure 25 
shows, for the £l — 1 model, the contribution of 
the different sed types to the background. The 
dominant contribution is from the cirrus compo- 
nent at most wavelengths, so the prediction is that 
more of the energy from starbursts is deposited in 
the general interstellar medium of a galaxy than is 
absorbed in the early stages close to the location 
of the massive stars. This dominance by the cir- 
rus component at submillimetre wavelengths also 
implies that many of the detected sources should 
turn out to be rather extended (kiloparsec scales 
rather than the more compact scales expected for 
nuclear starbursts). Table 5 gives the contribu- 
tions of the four components, and totals, to the 
background flux in the wavelength ranges 0.1-5, 
5-100 and 100-1250 /im. The totals are somewhat 
lower than those given by Gispert et al (2000), 
based on summing all claimed measurements of 
the background. 

7. Effect of including density evolution 

In the framework of hierarchical, bottom-up 
galaxy formation scenarios, like those based on 
cold dark matter, we expect that galaxies form as 
the result of mergers between smaller fragments. 
Thus we might expect to see a higher total density 
of galaxies as we look back towards the past. To 
test whether attributing part of the star formation 
history to density evolution has a major effect on 
the predicted counts and backgrounds, I have con- 
sidered, for the Q Q = 1 simple modification 
in which the comoving density of galaxies varies 
with rcdshift as 

p(z) = p(0)(l + z) n . 

For n=l, this means that the comoving number- 
density of galaxies at redshifts 2, 5 and 10 (our 
assumed cutoff redshift), is increased by a factor 
3, 6 and 11 compared with the present. This would 
represent a very substantial degree of merging of 
galaxies over the observable range of redshifts. 



Since the background intensity depends on the 
product of the luminosity and density evolution 
rates, the background spectrum will be unaltered 
if we simply increase P by 2n/3. At low redshifts 
the counts will be unaffected if P is increased by 
n/3. I find that for a combined density evolution 
model with n = 1 and luminosity evolution with P 
= 1.5, Q = 5.4 (and Sl = 1) the fits to the counts 
at 15-850 p,m hardly change over the observed 
range, but the predicted background is raised by 
about 0.2 dex, giving better agreement with ob- 
servations for this cosmological model. For the 
A=0.7 model the fit to the background would be 
significantly worsened for models consistent with 
the observed counts. 

8. Predicted confusion-limited performance 
of infrared and submillimetre tele- 
scopes 

The source-count models derived here can be 
used to compare the ultimate performance of dif- 
ferent infrared and submillimetre telescopes. To 
do this in a simple and transparent way, I treat 
confusion in a simplified but consistent manner. 
I assume that the 5-a confusion limit falls at a 
source-density of 1 source per 40 beams, defined 
to be circular apertures of radius 

Oconf = 0.6X/D 

where D is the diameter of the telescope. 
The corresponding source-density limit is denoted 
Nconfi P er sc l deg. The corresponding confusion 
limit, and the fraction of sources in different red- 
shift bins can then be estimated. The above 
assumption is valid for a population of sources 
with N(S) oc S' 1 - 5 (Murdoch et al 1973, Condon 
1974). The exact calculation depends slightly on 
the source-count slope. 

Table 6 gives the wavelength, 2 8 con f, N con f, 
Sconf, and redshift distribution for my best o = 
1 model, with (P,Q) = (1.2,5.4). 

9. Discussion and Conclusions 

(1)1 have developed a parameterized approach 
to the star formation history, which is sufficiently 
versatile to reproduce many proposed model histo- 
ries. The model assumes that the evolution of the 
star formation rate manifests itself as pure lumi- 
nosity evolution. I have stressed the importance 
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of ensuring that the assumed luminosity function 
is consistent with available 60 ^m rcdshift sur- 
vey data and that the assumed spectral energy 
distributions are realistic. The observed far ir 
and submm counts and background then provide 
strong constraints on the model parameters. The 
best fit for the fl = 1 model is close to the star- 
formation history predicted by the CDM scenario 
of Cole et al (1994), but is significantly less peaked 
than the predictions of Pei and Fall (1995). 

The models consistent with infrared and sub- 
millimetre counts and backgrounds tend to show 
a flat star formation rate from z = 1-3, consis- 
tent with the observed star formation history de- 
rived from HDF galaxies, using photometric red- 
shifts, and other uv, infrared and submillimetre 
data (these will be discussed in a subsequent pa- 
per, Rowan-Robinson 2000). The most striking 
difference from previous modelling work is the 
dominant role predicted for the cirrus component 
at submillimetre wavelengths. This is a conse- 
quence of the requirement that all galaxies share 
a common star formation history. The models 
presented here are by no means unique. If a 
very strongly evolving, heavily obscured popula- 
tion of starbursts is postulated, as in the models 
of Franceschini et al (1997) and Guiderdoni et al 
(1998), then such a population will affect only the 
submillimetre counts and background (and per- 
haps the faintest 60 /xm counts and redshift distri- 
butions). Existing 850 /im surveys (Hughes et al 
1998, Barger et al 1998, Eales et al 1999) suggest 
that at least some of the submm galaxies may be 
luminous, opaque starbursts. The prediction here, 
however, is that most will show a cooler sed, and 
greater spatial extent, than expected for a typical 
M82- or A220-like starburst. 

Areas requiring further work are (i) the need 
to consider the evolution of the shape of the seds, 
particularly for the cirrus component, with red- 
shift. The increased star- formation rate at earlier 
times would tend to make the dust temperature 
higher, but this is partially offset by the much 
lower abundance of low-mass stars at earlier times. 
The evolution of the metallicity may also affect 
the grain properties and hence the seds. (ii) the 
AGN dust tori models require a further parame- 
ter, the orientation, and this may affect the AGN 
counts at optical wavelengths, (iii) the tendency 
for the probability of finding an AGN component 



to increase with far infrared luminosity is not fully 
reflected in the approach followed here, (iv) it will 
be worthwhile to extend the predictions of these 
models to radio and X-ray wavelengths. 



9 



A. Cosmological formulae 

Collecting together the expressions for the cosmological observables, most of which are analytic functions, 
for the models consider here, we have, for (a) A = : 

the comoving radius 

r(z) = 2 (i - n )V2[n o2 + (n - 2)((i + n «)V2 - i)]/n 2 (i + z ), 

the luminosity distance 

D lum (z) = (1 - O,) 1 ^ (i +z ) r ( z ) ( c /tf o ), 

the comoving volume 

V(z) = 0.5 (r(l + r 2 ) 1 ^ + m (r + (1 + r 2 ) 1 / 2 ) (c/H ) 3 , 
and the age of the universe (in units of the Hubble time) 

H t(z) = [(1 - n o )(l + ^)] 1 /2/(Q (l + z )) - 0.5 ln [(2 - ft + Q z + 2[(1 - o )(l + fi z)] 1 /2)/(fi (l + *))], 

i/ t = (i - a,) 1 / 2 /^, - 0.5 in[(2 - o + 2(1 - n )V2)/n ]. 

For (b), k = 0, all the relations are analytic except that for r(z), which is an elliptic integral and has to 
be evaluated numerically. 

the comoving radius 

r(z) = 

2(1 - no) 1 ' 2 f [n /(i - n„) + y e r 1/2 dy 

the luminosity distance 
D lum (z) = (1+z) r(z) (c/# D ), 
the comoving volume 
V(z) = r 3 (c/H ) 3 /3, 

and the age of the universe (in units of the Hubble time) 
H t(z) = 2sinh- 1 [((l - r! o )/0 o ) 1 /2(i + z)- 3 / 2 ]/[3(l - a,) 1 /?], 
H t = 2sinh- 1 [{{\ - o )/r! o ) 1 /2]/[3(i _ Q )V 2 }. 
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Fig. 1. — Examples of star formation histories of 
the form eqn (1): (a) Fitted to those predicted by 
Pei and Fall (1995) (from the top, (P,Q) = (2.83, 
8.15), (3.23, 10.25), (3.52, 11.92)) from an infall 
model (their Figs Id, 3d, 4d), using dereddened 
quasar Lyman a absorber data, and (dotted curve) 
the prediction by Cole et al (1994) derived from a 
CDM scenario, (P,Q) = (1.5, 5.). (b) For n o = 1 
models consistent with Lilly et al (1996) CFRS 
data, defined to have </>*(l)/</>* (0) = 2 39 . From 
top at high z: (P,Q) = (1, 5.8), (2, 7.4), (3, 9.0), 
(5, 12.2). (c) Best-fit models for far infrared and 
submm counts and background for o = 1, (solid 
curve, (P,Q) = (1.2,5.4)); fl a = 0.3, (dotted curve, 
(P,Q) = 2.1, 7.3)); A = 0.7, (broken curve, (P,Q) 
= (3.0, 9.0)). The curves for (b) and (c) have been 
displaced upwards by +1 and +2 respectively. 

Fig. 2.— S(100)/S(60) versus 60 fim luminosity 
for PSCz galaxies. The solid line shows the trend 
with luminosity of the weighted average of the four 
component seds adopted. 



et al 1998), and log N850(l mJy) = 3.90 ±0.15 
(Blain et al 1999b). The locus corresponding to lg 
N60(50 mJy) = 1.30 (Hacking and Houck 1987) is 
almost identical to the 60 /jm locus, with a slightly 
wider uncertainty (± 0.28 in Q). 



Fig. 3.— S(25)/S(12) and S(60)/S25) versus 60 
fim. luminosity for PSCz galaxies. 

Fig. 4.— S(60)/S(850) versus 60 luminosity 
for galaxies observed by Rigopoulou et al (1996), 
Dunne et al (2000, open circles) and Fox et al 
(2000). 

Fig. 5. — Adopted spectral energy distributions 
for the four components adopted in this study: cir- 
rus (with optical emission split into low-mass (bro- 
ken curve) and high-mass (dotted curve) stars), 
M82-starburst, A220-starburst (models from Efs- 
tathiou et al 2000), AGN dust torus (model from 
Rowan- Robinson 1995), showing assumed opti- 
cal/ir ratio at log L60 = 14 (upper curve) and 
8. 

Fig. 6. — Average sed as a function of 60 fim lu- 
minosity, ranging from logi^L^/ '£©) = 8 to 14. 
The relative proportions of the different compo- 
nents, at 60 fim, are given in Table 1 and the av- 
erage is calculated using eqn (4) (AGN dust torus 
component excluded). 

Fig. 7. — P-Q diagram, with loci for models fitting 
60 /zm counts and V/V m test (solid curves) and 
850 /im counts (dotted curves). Counts loci corre- 
spond to log N850(4 mJy) = 2.95 ±0.15 (Hughes 
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Fig. 8. — Luminosity functions at 60 /xni for 
the 4 spectral components. Units of are 
(Mp^^dex^ 1 , luminosity (vL v ) in solar units. 
All luminosity functions arc for il — 1 model, H Q 
= 100 kms~ 1 Mpc~ 1 . Observed points are derived 
from PSCz data. 

Fig. 9. — Luminosity functions at 850 /im for the 
4 spectral components. The data are from Dunne 
et al (2000). 

Fig. 10. — Luminosity functions at 12 /im for the 4 
spectral components. Observed points taken from 
Rush et al (1993) (filled circles: Scyfcrts, open 
triangles: non-Seyferts). 

Fig. 11. — Luminosity functions at 0.44 /im for the 
4 spectral components. Data for quasars derived 
from PG sample and for galaxies from Loveday et 
al (1992). 

Fig. 12. — Bolomctric luminosity functions for 
the 4 spectral components. Note that for lumi- 
nosities greater than 10 12 5 L Q there is an approx- 
imately equal contribution from M82-type star- 
burst, Arp220-like starbursts and AGN dust tori. 

Fig. 13. — Integral source counts at 850 /im. Data 
are from Hughes et al (1998), Eales et al (1998), 
Smail et al (1997), Barger et al (1999), Blain et 
al (1999c), Fox (2000). The 3 models shown are, 
from bottom at faint fluxes, for Q a = 1 and (P,Q) 
= (1.2, 5.4), (solid curve), for Q Q = 0.3, (P,Q) = 
(2.1, 7.3) (dotted curve) and for A = 0.7, (P,Q) = 
(3.0, 9.0) (broken curve). 

Fig. 14. — Source counts at 175 /im. Data points 
are from Kawara et al (1998), Guiderdoni et al 
(1998), Dole et al (2000). Models as in Fig 13. 
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Fig. 15. — Source-counts at 90 /im. Data from 
IRAS PSCz (triangles) and ELAIS (filled circles) 
are from Efstathiou et al (2000). Models as in Fig 
13. 

Fig. 16. — Source counts at 60 /im. Data arc from 
Lonsdale et al (1990) (at 0.2-10 Jy), Hacking and 
Houck (1987) (at 50-100 mJy), Rowan- Robinson 
et al (1991), Gregorich et al (1995) (higher point 
at 50 mJy: see Bertin et al 1997 for discussion of 
possible contribution of cirrus to Gregorich et al 
counts). Models as in Fig 13. 

Fig. 17. — Source counts at 15 /urn. Data 
from Oliver et al (1997, O), Serjeant et al (2000, 
ELAIS), Elbaz et al et al (2000, E) Rush et al 
(1993, R, crosses), Verma (2000, V, filled circles). 
Models as in Fig 13. The lower dash-dotted line 
shows counts of the AGN dust torus population 
for the Cl a = 1 model. 

Fig. 18. — Differential source counts at 2.2 /im. 
Data from McCracken et al (2000). Models as in 
Fig 13. 

Fig. 19. — Differential source counts at 0.44 /im. 
Galaxy data is from Metcalfe et al (1995), quasar 
data from Boyle et al (1988). Models as in Fig 13: 
Lowe dotted curve is prediction for quasars in Q Q 
= 1 model. The long-dashed line shows the effects 
of including density evolution (see text) on the il 
= 1 model. 

Fig. 20. — Predicted redshift distribution at 850 
/im, logioS(Jy) = —2.7. Bin centred at z = 5.25 
refers to z > 5. 

Fig. 21. — Predicted redshift distribution at 175 
F, log w S(Jy) = -1.0. 
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Fig. 22. — Predicted redshift distribution at 60 
/zm, logwS(Jy) = —0.7. Contribution from cirrus 
is shown dotted, that from starbursts shown as 
broken histogram. Observed points are derived 
from the FSSz redshift survey of Oliver et al 1996. 

Fig. 23. — Predicted redshift distribution at 0.44 
/im, logiaS(Jy) = —7.0 (corresponding to Bab 
= 26.4). Dotted and dash-dotted show contri- 
butions of high-mass and low-mass stars, broken 
curve shows contribution of starburst component. 

Fig. 24. — Predicted spectrum of integrated back- 
ground for same models as Fig 9. Data from 
Fixsen ct al (1998) (far ir and submm), Pozzetti 
ct al (1998) (opt and uv), Serjeant et al (2000) (15 

Fig. 25. — Predicted spectrum of integrated back- 
ground for o = 1 model, showing contribution of 
the different components. The contribution of the 
Arp220-like starbursts is less than 0.01 nWm~ 2 
at all wavelengths. The upper solid curve at long 
wavelengths shows the effect of including density 
evolution (see text). 
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Table 1: ti(Lea), proportion of different sed types 
at 60 (im, as function of Lqq. 



IqioLso/Lq 


cirrus 


M82-starburst 


A220-starburst 


AGN dust torus 


8.0 


0.99 


0.001 


0.5e-10 


0.0091 


9.0 


0.80 


0.1988 


0.5e-7 


0.0012 


10.0 


0.50 


0.4999 


0.5e-4 


0.5e-4 


11.0 


0.10 


0.85 


0.05 


0.19e-4 


12.0 


0.1e-3 


0.75 


0.25 


0.25c-6 


13.0 


0.1e-7 


0.75 


0.25 


0.25e-8 


14.0 


O.lc-12 


0.75 


0.25 


0.25e-10 



Table 2: Observed source-counts 



wavelength (^m) 


logio S(Jy) 


log 10 N(> S) 


reference 


15 


-3.60 


3.56 ± 0.135 


Oliver et al (1997) 




-2.92 


1.93 ± 0.13 


Elbaz et al (2000) 




-3.30 


2.93 ± 0.13 


Elbaz et al (2000) 




-3.70 


3.63 ± 0.13 


Elbaz et al (2000) 




-4.00 


3.90 ± 0.09 


Elbaz et al (2000) 


60 


1.00 


-2.12 


Lonsdale et al (1990) 




0.5 


-1.39 


Lonsdale et al (1990) 




0.1 


-0.85 


Lonsdale et al (1990) 




-0.22 


-0.38 


Rowan- Robinson et al (1991) 




-0.60 


0.21 ± 0.1 


Lonsdale et al (1990) 




-0.70 


0.38 


Lonsdale et al (1990) 




-1.00 


0.76 ± 0.07 


Hacking and Houck (1987) 




-1.30 


1.46 


Gregorich et al (1995) 




-1.30 


1.30 ± 0.08 


Hacking and Houck (1987) 


175 


-1.0 


1.59 ± 0.1 


Kawara et al (1998),Guiderdoni et al (1998) 


850 


-2.40 


2.95 +0.15,-0.25 


Hughes et al (1998) 




-2.70 


3.40 +0.15,-0.25 


Hughes et al (1998) 




-2.55 


3.25 +0.13,-0.17 


Eales et al (1998) 




-2.40 


3.39 ± 0.19 


Smail et al (1997) 




-2.65 


3.37 ± 0.30 


Barger et al (1999) 




-3.00 


3.90 +0.14,-0.21 


Blain et al (1999c) 




-2.155 


2.56 +0.23,-0.52 


Fox (2000) 
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Table 3: Observed submillimetre background 



wavelength (//m) 


log\ovl v (nW/ m 2 / sr) 




reference 


140 


1.51 +0.14, -0.23 


Schlegel et 


al 


(1998) 




1.40 +0.11, -0.15 


Hauser et 


al 


(1998) 




1.10 ± 0.10 


Fixsen et 


al 


(1998) 




1.18 +0.15, -0.23 


Lagache et 


al 


(1999) 


240 


1.00 ± 0.20 


Puget et 


al 


^ryyo ) 




1.23 +0.09, -0.12 


Schlegel et 


al 


(1998) 




1.15 +0.08, -0.11 


Hauser et 


al 


(1998) 




1.04 ± 0.10 


Fixsen et 


al 


(1998) 




1.06 +0.06, -0.08 


Lagache et 


al 


(1999) 


350 


0.85 ± 0.20 


Puget et 


al 


(1996) 




0.78 ± 0.10 


Fixsen et 


al 


(1998) 




0.85 +0.13, -0.21 


Lagache et 


al 


(1999) 


500 


0.37 


Fixsen et 


al 


(1998) 




0.59 +0.16, -0.26 


Lagache et 


al 


(1999) 


750 


-0.12 


Fixsen et 


al 


(1998) 




-0.05 +0.11, -0.14 


Lagache et 


al 


(1999) 


850 


0.00 ± 0.15 


Puget et 


al 


(1996) 




-0.30 ± 0.10 


Fixsen et 


al 


(1998) 



Table 4: Best fit star-formation history models for 
n o = I/O, = 0.3/A = 0.7 



cosmological P Q igo £ ^* 

model hL oc [otMpc~ 3 



n Q = 1 


1.2 


5.4 


4.3 


5.70 


0.0027 hr 


i 


n D = 0.3 


2.1 


7.3 


4.4 


6.66 


0.0032 hr 


l 


A = 0.7 


3.0 


9.0 


4.1 


7.25 


0.0033 hr 


l 



Table 5: Predicted integrated background flux 
from different components, for f2 = l/f2„ = 
0.3/A = 0.7 



component 



cirrus+star light 
M82 starburst 
A220 starburst 
AGN + dust torus 



vl v at 0.1-5 fim 

(nW/m 2 /sr) 



13.81/15.73/18.07 
1.25/1.71/2.11 
0.0004/0.0006/0.0007 
0.243/0.320/0.432 



5-100 urn 



5.30/6.94/9.01 
2.68/3.56/4.35 
0.004/0.005/0.006 
0.247/0.333/0.431 



100-1600 urn 



10.09/14.30/18.45 
1.83/2.63/3.25 
0.0074/0.0103/0.0122 
0.0074/0.0106/0.0116 



total 



28.20/36.97/45.36 
5.76/7.90/9.71 
0.0118/0.0159/0.0189 
0.490/0.674/0.875 



total(all types) 14.34/17.75/20.61 8.21/10.86/13.63 11.91/16.95/21.73 34.46/45.56/55.97 
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Table 6: Confusion limit for different infrared and 
submillimetre telescopes, for Sl Q =1 model 



Telescope D(m) X(fim) 2 9 conf N conf lg S conf 

(/sq dg) (Jy) 





0.0 


1 c 

15 


6.2 


10800 


A 11 

-4.zz 






90 


37.1" 


300 


-1.66 






175 


72.2" 


79 


-0.97 


SCUBA 


15. 


450 


7.4" 


590 


-1.90 






850 


14.0" 


2100 


-2.63 


SIRTF 


0.85 


8.0 


2.33" 


75780 


-6.17 






25 


7.3" 


7760 


-3.87 






70 


20.4" 


990 


-2.33 






160 


46.6" 


190 


-1.23 


ASTRO-F 


0.6 


60 


25" 


700 


-2.36 






175 


72.2" 


79 


-0.97 


FIRST 


3.5 


90 


6.4" 


10185 


-3.13 






120 


8.5" 


5730 


-2.50 






175 


12.4" 


2680 


-1.96 






250 


17.7" 


1320 


-1.73 






350 


24.8" 


670 


-1.70 






500 


35.4" 


330 


-1.78 






600 


42.5" 


230 


-1.88 


PLANCK 


(1.5) 


350 


5' 


4.6 


-0.83 






550 


5' 


4.6 


-1.13 






850 


5' 


4.6 


-1.63 



rcdshift distribn (fraction) 





n i 

0-1 




1 1 

1-2 




2-3 


3-4 


4-5 


>5 





,7011 


0, 


.2858 





.0077 


0.0009 


0.0001 


0.0 


0. 


.9708 


0. 


,0290 





.0002 


0.0 


0.0 


0.0 


0. 


ATO C 

.9785 


0. 


0212 





Anno 

,0002 


A A 

0.0 


A A 

0.0 


A A 

0.0 


0. 


.3287 


0. 


.4850 





.1583 


0.0258 


0.0020 


0.0002 


0. 


1082 





,3193 





.2801 


0.1772 


0.0746 


0.0406 





.3829 


0, 


,4525 





.0868 


0.0413 


0.0172 


0.0191 


0. 


,6062 


0. 


.3484 





.0428 


0.0024 


0.0001 


0.0 


0. 


,9314 


0, 


.0675 





.0010 


0.0 


0.0 


0.0 


0. 


,9741 


0. 


.0253 





.0003 


0.0 


0.0 


0.0 


0. 


.9334 





,0655 





.0010 


0.0 


0.0 


0.0 


0. 


.9785 


0. 


,0212 





.0002 


0.0 


0.0 


0.0 


0. 


.7456 


0, 


,2316 





.0211 


0.0016 


0.0001 


0.0 


0. 


,8366 


0, 


,1546 





.0083 


0.0004 


0.0 


0.0 


0. 


.8261 


0, 


,1683 





.0052 


0.0003 


0.0 


0.0 





.7302 


0. 


,2609 





.0085 


0.0004 


0.0 


0.0 


0. 


.5469 


0. 


.4077 





.0432 


0.0021 


0.0001 


0.0 





,3454 


0. 


,4838 





.1469 


0.0223 


0.0016 


0.0001 


0. 


.2426 


0. 


.4805 





.2145 


0.0544 


0.0071 


0.0010 


0. 


,8104 


0. 


,1851 





,0043 


0.0001 


0.0 


0.0 


0. 


.4454 


0, 


.4220 





.1222 


0.0097 


0.0005 


0.0001 


0. 


,1555 


0, 


,3732 





.3086 


0.1310 


0.0276 


0.0040 
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